AD  605958 


iltG.  NO  -— 
I.OG.NO - 


yjDSn. 


THE  RAMO-WOOLDRIDGE  CORPORATION 
Qxuded  Ivlissile  Research  Division 


I 

^'5 


ON  THE  SLOSHING  OF  LIQUID  ^ 
IN  A  CYLINDRICAL  TANK  J 

by 

John  W.  Miles 

^  J 

Report  No.  AM  6-5 


COPY  _ _ L 

HARD  COPY 

1 .  A  cn^ 

MICROFICHE 

%.  C  <S'o 

20  April  1956 


Aeromechanics  Section 


G1,*-TR-18 


M,  V.  Barton 


SUMMARY 


The  free  oscillations  of  a  slightly  viscous  liquid  in  a  cylindrical 
tank  (of  arbitrary  cross  section)  with  a  free  surface  are  analyzed  and 
the  decrement  of  motion  caused  by  laminar  boundary-layer  friction 
and  time-varying  depth  calculated.  An  increasing  or  decreasing  (in 
time)  depth  is  shown  to  have  a  destabilizing  or  stabilizing  effect, 
respec  .vely,  but  the  magnitude  of  this  effect  is  negligible  for 
practical  configurations.  The  results  are  applied  to  a  circvilar  tank, 
and  it  is  shown  that  the  depth  may  be  considered  infinite  when  it 
exceeds  the  diameter.  The  sloshing  oscillations  in  a  10-foot 
(diameter)  tank  are  found  to  have  a  (longest)  period  of  1.  8  seconds 
and,  assuming  the  kinematic  viscosity  to  be  that  of  water  at  20° C, 
a  time  to  damp  to  one -half  amplitude  of  7  minutes  for  mean  amplitudes 
smaller  than  about  10  inches.  The  transverse,  oscillatory  force  on 
the  tank  walls  associated  with  a  sloshing  amplitude  of  one  foot  would 
have  an  amplitude  of  2700  pounds  and  would  act  2,  7  feet  below  the 
mean  surface. 


SECTION  I  -  INTRODUCTION 


We  shall  consider  the  free  oscillations  of  a  slightly  viscous 
fluid  having  a  free  surface  but  otherwise  constiained by  a  vertical, 
cylindrical  tank  of  varying  (in  time)  depth.  The  preliminary  analysis 
will  be  carried  out  for  a  cylinder  of  arbitrary  cross  section,  after 
which  the  results  will  be  applied  to  fuel  sloshing  in  a  circular  tank.  * 
The  effects  of  friction  on  the  motion  of  the  liquid  are,  by 
hypothesis,  small  and  are  confined  to  a  thin  boundary  layer  at  the 
wall.  It  follows  that  a  first  approximation  to  the  flxiid  motion,  adequate 
except  for  the  calculation  of  damping,  may  be  calctdated  by  neglecting 
friction  entirely;  this  approximation  then  may'be  used  as  a  basis  for 
the  calculation  of  the  boimdary  layer  flow.  We  find  it  expedient  in 
this  approach  to  aim  at  an  energy  formulation  and  toward  this  end 
consider  the  kinematics  of  this  first  approximation  in  Section  II, 
establish  the  corresponding  kinetic  and  potential  energies  in  Section  III, 
find  the  laminar  boundary  layer  solution  and  corresponding  dissipation 
fiinction  in  Section  IV,  and  calculate  the  frequency  and  decrement  of 
the  mean  motion  in  Section  V.  The  local  equation  of  motion  does  not 

I 

appear  throughout  this  sequence  but  is  introduced  a  posteriori  in 
Section  VI  to  calculate  the  dynamic  forces  on  the  tank  walls.  he 
resvilts  established  in  Sections  II  through  VI  are  applied  to  rcular 
tank  in  Section  VII  and  a  numerical  example  considered. 


This  problem  also  has  been  considered,  from  a  somewhat  different 
viewpoint,  byR.R.  Berlot  in  Ramo- Wooldridge  Report  GM- TN-2 
(22  March  1956). 
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SECTION  n  -  KINEMATIC  FORMULATION 


We  consider  a  small  disturbance,  ^  (x,  y,  t),  of  a  free  surface  of 
liquid  from  its  eqmlibrium  position  z  =  0  in  a  cylindrical  tank  of  cross 
section  A  and  depth  h,  as  shown  in  Figure  1.  The  effects  of  viscosity 
being  negligible  in  first  approximation,  the  fluid  motion  will  be 
irrotational,  in  consequence  of  which  we  may  express  the  velocity  at 
any  point  in  the  liquid  as  the  negative  gradient  of  a  velocity  potential  4». 

The  compressibility  of  the  liquid  also  being  negligible,  this  potential 
must  satisfy  Laplace's  equation 

V^4,  =0  (2.1) 

The  kinematic  boundary  conditions  on  <}),  as  dictated  by  the  assxuned 
velocity  at  the  free  surface  and  the  condition  of  zero  normal  velocity 
at  the  walls  of  the  tank,  are 


on 

z  =  0 

\Z,  2a) 

0 

on 

z  =  -h 

(2.2b) 

0 

on 

S 

(2.2c) 

where  (2.2a)  is  imposed  at  the  mean  position  of  the  free  surface  by 
virtue  of  the  assumption  of  small  displacement;  S  denotes  the  lateral 
area  of  the  cylinder,  and  n  is  the  inwardly  directed  normal  to  S. 


Figure  1.  Cross  Section  of  Cylindrical  Tank 
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The  solution  to  (Z,  1)  and  (2.2),  due  originally  to  Poisson  and 

2  3 

Rayleigh,  may  be  exhibited  in  the  form 

i  (x,  y,  t)  =  q(t)  »|j(x,  y) 

<jj  (x,  y,  z,  t)  =  -qCt)  4'(x,  y)  cosh[[k(z+h)]  /k  sirih  (l<h:)  .  | 


(2.3) 

(2.4) 


where  4*  is  a  solution  to  the  Helmholtz  equation 

y  ^  4i+k^4j  =  0  (2.5) 

subject  to  the  boundary  condition 


i}^  =  0  on  C  .  (2.6) 

H 

where  C  denotes  the  boundary  of  the  cylinder  in  a  plane  of  xonstant  z. 

The  homogenous  boundary-value  problem  presented  by  (2i-5} 

and  (2.6)  is  analogous  to  that  presented  by  a  vibrating  membrane  or 

4 

propagation  of  sound  in  a  Cylindrical  tube  and  possesses  a  doubly 
infinite  number  of  solutions,  say  4^^..,  corresponding  to  those  values 
of  the  parameter  k,  say  k^j,  that  are  permitted  by  the  boundary  condition 
(2.  6),  The  4^  as  is  well  known,  are  orthogonal,  so  that  the  correspond- 
ing  modes  of  motion  are  independent;  accordingly,  we  may  omit  the 
subscripts  with  the  implication  that  4^  denotes  one  mode  only.  We 
find  it  convenient,  in  the  subsequent  analysis,  to  render  4^  dimension¬ 
less  and  normalize  the  individual  modes  such  that 


dA  =  A 


(2.7) 


We  also  have,  from  the  consideration  (dictated  by  continuity)  that  the 
mean.,  elevation  represented  by  ^  must  vanish,  the  relation 


dA  =  0 


(2.8) 


We  emphasize  in  concluding  this  section  that  the  formulation  to 
this  point  rests  on  purely  kinematical  assumptions,  viz. ,  irrotationality 
and  .incompre  s  sibility . 


3 


SECTION  m  -  KINETIC  AND  POTENTIAL  ENERGIES 


The  function  q(t)y  dimensionally  a  length,  may  be  regarded  as  a 
generalized  coordinate  for  the  particialar  mode  {ij^and  k)  to  which  it 
corresponds.  The  kinetic  energy  of  this  mode,  allowing  for  variations 
of  depth  with  time,  is  given  by 


(3.1a) 


(3.1b) 


where  p  denotes  the  density  of  the  liquid  and  m  the  mass  per  imit 
depth  (i.e.,  m  =  pA).  The  individual  contributions  of  the  three  terms 
in  the  integrand  of  (3.  lb)  may  be  evaluated  according  to 


m 

ZXJ 


n 


rc  •  2 

j  A  ^  j  mh 

dA  h  dz  = 

J 


(h  +  ^)  dA  =  i  mhh^ 


(3.2) 


m 

2X 


dA 


u 


-h 


(-2h4>^)  dz  = 


=  X'  M«t‘^  =  n'hqq 


z=“h 

z=0  ^ 

z=0 


dA 


(v^) 


ri 


z=0 


4>4>2  +  0(r) 


dA  = 


coth  (kh) 


.2 

q 


(3.3) 

(3.4) 


where  those  integrals  of  0{t,)  vanish  identically  in  virtue  of  (2.8)  and 
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those  of  0(^  )  have  been  evaluated  with  the  aid  of  (2.  7/;  the  volume 
integral  in  (3.4)  has  been  converted  to  a  surface  integral  with  the  aid 
of  Green's  theorem,  (2,  1),  and  (2.2).  The  end  result  for  the  kinetic 
energy  then  reads 


T  - 


T 

o 


•  i 

mhqq  +  jm 

1  '  2 
Y  (mh)  h 


coth  .(kh)  ^  •  2 


(3.5a) 
(3.  5b) 
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n 


where  is  the  kinetic  energy  associated  solely  with  the  changing  level 
of  the  liqviid,  while  T  —  represents  the  additional  contribution  of  one 
mode  only;  more  generally,  the  terms  on  the  right-hand  side  of  (3.  5a) 
may  be  siimmed  over  all  modes  but  do  not  involve  cross  products  of 
the  generalized  coordinates. 

The  potential  energy  relative  to  the  bottom  of  the  tank  is  given 


by 


--  {3.6a) 


(h^  +  2h;  +  dA 


(3.6b) 


Evaluating  the  integrals  with  the  aid  of  (2.7)  and  (2.8),  we  obtain 


V  -  =  i  mg  (3.  7a) 

=  i  mg  (3.  7b) 

where  denotes  the  potential  energy  of  the  imdisturbed  fluid.  We 
remark  that  g,  the  acceleration  of  gravity,  shovild  be  replaced  by  an 
equivalent  value  (say  ng)  if  the  entire  tank  is  subjected  to  an  acceleration. 

We  shall  proceed  further  on  the  assumption  that  h  varies  slowly. 

The  motion  q(t)  then  will  be  approximately  simple  harmonic,  arid  the 
mean  values  of  the  perturbation  energies  may  be  approximated  by 
(neglecting  the  variation  of  h  in  averaging  over  one  cycle) 


1  coth  (kh) 
k 


.  2 
q 


(3.8) 


(3.9) 


SECTION  TV  -  VISCOUS  DISSIPATION 


We  shall  calculate  the  frictional  dissipation  of  the  free  oscillations 

on  the  assumption  that  viscous  forces  are  negligible  everywhere  except 

in  a  laminar  boundary  layer  at  the  walls  of  the  tank,  reflecting  the 
5 

well-known  fart  that  the  logrithmic  decrements  associated  with  friction 
inside  and  outside  of  this  layer  are  0(v/cjll  )  and  0(v/a)R  )  , 
respectively,  where  v  is  the  kinematic  viscosity,  co  the  angular  frequency, 
and  R  a  characteristic  length  proportional  to  area/perimeter  (often 
called  the  hydraulic  radius  ) .  The  thickness  of  the  boundary  layer  being 
small  compared  with  all  other  characteristic  lengths,  we  may  also 
treat  the  viscous  flow  as  plane  and  neglect  the  pressure  gradient.  The 
linearized  equation  for  the  velocity  v  in  the  bovuidary  layer  then 
contains  only  the  inertial  and  shear  terms  and  reads ^ 


(4.1) 


We  reqixire  a  solution  to  this  eqviation  that  vanishes  outside  the  boundary 

layer  (n — ►  oo)  and  just  cancels  the  tangential  velocity,  say 

v^^  cos  {(Jit  +  €),  calculated  at  the  wall  in  the  absence  of  friction,  viz. , 

~v - ►-  0  ,  n - ►-  oo  (4.  2a) 

v^  =  cos  (wt  +  c)  ,  n  =  0  (4.2b) 

The  problem  posed  by  (4. 1)  and  (4.2)  corresponds  to  that  for  an 

7  8 

oscillating  flat  plate,  first  solved  by  Stokes  with  the  result 

=  e  cos  (u)t  -  pn  +  €)  (4.  3a) 

p=(a)/2v)^^^  (4.3b) 
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The  integrate  !  flow  defect  per  unit  width  of  the  boundary  layer  is 


00 


J  o 


dn  = 


V  C08  (u)t  -  Tr/4) 
o 


'ITp 


Rj  = 


6 

o  1  !  _ 

1  o 

A 


where 


(4.4) 


We  infer  from  this  result  that  the  mean  flow  defect  in  the  boundary 
layer  lags  the  outer  flow  by  45^,  and  that  the  displacement  thickness 

-  6,  =  — ! 14.  5) 

‘  'I  “ 

The  aorresponding  Reynolds  number  is 


\JZ  V  vfZuTv 

and  A(=«”*  j  denote  the  r.m.  s.  velocity  and 

amplitude  of  the  motion  outside  the  boundary  layer. 

The  question  of  the  stability  of  the  laminar  flow  given  by  (4.  3) 
does  not  appear  to  have  been  examined  in  the  literature,  but  we  may 
approximate  the  amplitude  of  motion  necessary  to  initiate  turbulence 
by  comparison  to  the  transition  problem  for  steady  flow  along  a  semi- 
infinite  flat  plate.  The  minimum  transition  Reynolds  number  for  the 

■Oi 

steady  flow  problem,  based  on  displacement  thickness,  is  found 
experimentally  to  vary  between  about  500  and  1700,  depending  on  the 
level  of  free  stream  turbtilence^^  (the  theoretical  minimum  value  is 
about  420^^).  Taking  the  lowest  experimental  value,  we  infer  that  the 
oscillatory  flow  in  the  boundary  layer  will  remain  laminar  if 

A  <  500  6j  =  500,J~^ 


(4.6) 


(4.  7) 


We  remark  that  this  transition  Reynolds  number  is  rather  insensitive 
to  the  pressure  gradient  in  the  outer  flow9,  from  which  we  surmise 
that. the  detciils  of  the  outer  flow  are  sufficiently  unimportant  to  warrant 
a  moderate  expectation  of  "'uccess  in  the  extrapolation  to  an 
oscillatory  flow. 
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This  resvdt  also  may  be  obtained,  at  least  as  to  order  of  magnitude, 
on  the  basis  of  a  transition  time  (equal  to  transition  distance  divided 
by  free  stream  velocity  for  the  steady  flow  problem)  Reynolds  number 
if  the  allowabJe  time  in  the  oscillatory  problem  is  equated  to  some 
fraction  of  the  period;  the  numerical  value  of  (4.  7)  is  obtained  if  this 
fraction  is  assvimed  to  correspond  to  a  phase  shift  of  approximately  45°. 

The  viscous  shear  in  the  boundary  layer  is  given  by  -p  (9  v/^n), 
so  that  the  rate  of  dissipation  per  \init  area  is,  in  terms  of  a  dissipation 
function 


9F 


00 


9  V 

fflT 


dn 


u 

We  reqxiire  only  the  mean  value  of  this  dissipation  function,  viz. , 
^  9  F 


(4.8) 


9S 


=  2 


2 


00 


,.-2pn 


c 


-1  2 


cos(u)t  “Pn  +  € )  -  sin({ot  -  Pn  +  6  ) 


o 


00 


-Zpn  , 

e  dn 


dn  (4. 9a) 
{4.9b) 


u 

2 


(4.9c) 


We  consider  first  the  integral  of 

2 


over  the  bottom,  where 


=  2  q^  (yi|r)^/k^  sh^  (kh) 

|z  =  -h 

the  factor  of  two  accounting  for  the  fact  that  the  mean  square  of  v 
is  one-half  the  square  of  its  amplitude  in  consequence  of  the  assumed 
harmonic  motion.  Neglecting  the  variation  of  depth  over  one  cycle, 
we  then  have  «  ^ 


(4.  10) 


z=  -h 


dA  = 


,  .2 

2q 


k^  sh^  (kh)  vJ 


k^  sh^  (kh)''' 


C7ijj)  dA 


dA 


=  2A^^  csch^  (kh) 


(4.11a) 

(4. 11b) 
(4.  11c) 
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where  the  integral  hae  been  reduced  with  the  aid  of  Green's  theorem 
and  (2.5)  through  (2.7).  Similarly,  the  Integral  over  the  lateral  wall 


_|Z_  <fd  f 

•h*  (kh)  J  J  -H  [ 


h)J 


[>.(J5+h)] 


+h  csch^(kh) 


] 


(4.  12a) 


.  h  c.ch^^  (kh)J  +2 


(4.  12b) 


Adding  (4.  12b)  to  (4.  11c)  and  multiplying  the  total  by  |iP/2,  the  total 
diesipation  per  cycle  may  be  placed  in  the  form 


ZP  *  2  a  m  k"^  coth  (kh) 


where 


iq  '  X  *  (f  Ij)  ] 


(4.13) 


(4.  14) 


(4. 15a) 


(4. 15b) 


We  may  designate  the  lengths  Rj  and  R2,  which  are  proportional  to 
wetted  area  divided  by  perimeter,  as  modified  hydraulic  radii. 


SECTION  V  -  FREQUENCY  AND  DECREMENT  OF  MEAN  MOTION 


The  mean  (over  one  cycle)  motion  of  the  approximately  harmonic 
motion  will  be  governed  by  the  equation 

^  [cm;)  +  (5.1) 

provided  that  the  change  over  one  cycle  ii  small.  The  mean  values  of 
the  perturbation  energies  are  given  by  (3.8)  and  (3.9).  while  the 
dissipation  function  is  given  by  (4.9);  noting  that 

4  =  q  (5. 2) 

for  an  approximately  harmonic  motion,  these  means  become 


mw^  k”^  coth  (kh)  ^ 

“7 

(5.3) 

mg  q 

(5.4) 

zr  =  2 

a  m«^k’^  coth  (kh)*!^ 

(5.5) 

Substituting  (5.  3)  through  (5.5)  in  (5. 1)  and  dividing  through  by  mg, 
we  obtain 


(5.6) 


We  may  obtain  an  adequate  approximation  to  the  frequency  in 
(5.  6)  by  noting  that  the  mean  values  of  kinetic  and  potential  energy 
would  be  exactly  equal  for  a  truly  harmonic  motion,  whence 

tJ"  =  kg  tanh  (kh)  (5.7) 
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in  agreement  with  the  classical  result  for  constant  depth.  The 
time  dependence  of  w,  as  given  by  (5.  7),  and  the  deviations  from  this 
value  associated  with  viscous  dissipation  are  both  second-order 
effects  in  (5.  6),  so  that  we  may  substitute  (5.  7)  directly  in  (5,  6)  to 
obtain 

T  —K 

■  +  |^2a  -kh  csch  {2kh)J  q  =0  (5.8) 


which  yields 


2 

q 


exp 


-  kh 


csch 


a 


(5.9) 


where  is  the  initial,  mean  square  value  of  q. 

The  damping  ratio  (for  the  r.m.  s.  amplitude)  implied  by  (5.9)  is 


-1 

7  =  0) 


1  » 

a  -  y  kh  csch  (2kh) 


Substiuting  a  from  (4.  8),  we  have 


f  kh  csch  (2  kh) 


* 

h 

2o)h 


(5.10) 


(5.11) 


We  note  that  time  varying  depth  has  a  destabilizing  or  stabilizing  effect 

as  h  is  positive  or  negative,  respectively.  The  magnitude  of  the 

effect  is,  however,  rather  small  for  practical  values  of  (h/ojh). 

In  most  applications  kh  will  be  sufficiently  large  (say  kh  >  2) 

-2kh 

to  justify  the  neglect  of  terms  of  order  e  ,  and  the  results  (5.  7)  and 
(5.  11)  simplify  to 


(5.  12) 
(5.13) 


In  this  approximation  the  effect  of  varying  depth  is  negligible. 
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SECTION  VI  -  PRESSURE  DISTRIBUTION 


The  pressure  distribution  (p)  on  the  lateral  walls  of  the  tank 
may  be  calcvilated  from  the  Bernovilli  equation 

P  -  Pq  =  P(4>^  -  gz) 


where  p^  denotes  the  pressure  above  the  free  surface  (the  effects 
of  friction  and  changing  depth  on  the  pres  sure- are  neglected).  The 
time  derivative  of  the  potential,  as  given  by  (2.  4) ,  is 


=  ••  “I  ch  Qk(z  +  h)  J  /  k  8h(kh)  > 
^  - 

=  co^  ^  -^ch  ^k(z  +  h)J/  k  8h(kh)  | 
=  g  ^  |^ch£k{z  +h)j/  ch{kh)  ^ 


where,  in  (6.  2c),  w  has  been  substituted  from  (5.  7).’  Substituting 
(6,2c)  in  (6.  1),  we  obtain 


P  -  Pq  =  Pg 


^ch  [k(z  +  hj| 
ch(kh) 


-  z 


We  remark  that,  as  required,  p  =  p  at  z  =  ^  (neglecting  terms  of 
2  ° 

order  t  );  indeed,  this  requirement  affords  the  usual  derivation  of 
(5.7). 

The  integrated  pressure  on  a  vertical  strip  of  the  tank  is  given 

by 


-h 


(p  -  p^)  dz  =  p  g 


I 


^  sh  Ijk(z  4  h)] 
k  ch(kh) 


-h 


+  I  (h^  - 


1 


=  ipgh^  +  pgC  k'^  tanh(kh)+0(  C^) 


(6.1) 

(6.2a) 
(6.2b) 
(6. 2c) 

(6.3) 


(6.  4a) 

(6.4b) 
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The  corresponding  moment  with  respect  to  the  bottom  of  the  tank  is  given 
by 


(P  -  Pq)(z+1i)  <iz  =  pgh^  +  pg^,  '|k”^  h  tanh  (kh) 

+  k‘^  [  sech  (kh)  -  l]  j  +  O  ( 


(6.  5a) 


=  ^  pgh^  +  pg  4  k~  ^  tanh  (kh) 


•  h  - 


1 


1 


tanh 


(kh/2)J 


(6.  5b) 


The  center  of  pressure  of  the  perturbation  (^)  component  is,  therefore, 
a  distance 


d  =  h  -  k  ^  tanh  (kh/2) 


(6.6) 


from  the  bottom  of  the  tank. 
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SECTION  VII  -  CIRCULAR  TANK 


We  consider  now  the  application  of  the  foregoing  results  to  the 

more  specific  configuration  of  a  circular  cylinder  of  radius  a.  The 

approximate  solutions  of  (2.  5)  and  (2,  6)  read,  in  the  cylindrical  polar 

14 

coordinates  r  and  0, 


i{j{r,  0)  =  C  ,  J  (kr)  cos  (s0)  ,  s  =  0,  1,  2,  .  ,  . 

Sic  s 

where  denotes  Bessel's  fxinction  of  order  s,  and  k  is  any  root  to 


the  transcendental  equation 


Jg  (ka)  =  0 


The  normalizing  coefficient  C^^^  is  determined  by  the  requirement 
(2.7),  viz.. 


1  = 


TT  a 


2-ir 


J^ikr)  cos^  (s  0)  rd rd 0 


r,  2  2  2 

I  ^  ^  -  s 

2k 


r-2-  1  V 


The  dominant  mode  of  oscillation  corresponds  to  the  smallest 
root  of  (7.  2)  and  is  specified  by 


s  =  1  ,  ka  =  0 .  586  ir  =  1 . 84 

The  corresponding  frequency,  as  determined  by  (5.7),  is  given  by 

CJ  1 


2T  =z7\/^g  ™ 

=  0.216  y  I  tanh  (1. 84  h/a) 


(7.1) 

(7.2) 

(7.3) 

(7.4) 

(7.5a) 

(7.5b) 
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The.  (reciprocal)  hydraulic  radii,  obtained  by  substituting  (7.  1) 
and  (7.  3)  in  (4. 15a,  b),  are  found  to  be 


1  _  2  +  s^ 


(7.6a) 


1  ^1^  2 
^  h"a 


(7.  6b) 


The  corresponding  deunping  ratio,  obtained  by  substituting  (7,5a)  and 
(7.  6a,  b)  in  (5.11)  and  neglecting  (h/wh),  is 


-1/4,,  ,-1/4,  .1/4,,,, 

g  (ka)  tanh  (kh) 


+  (ka  -  2kh)  csch  (2kh) 


(7.7) 


which  reduces  further  to 

•  tanh^/^  (1.84h/a)  (7.8) 

for  the  dominant  mode  of  (7.4). 

The  depth  correction  functions  for  frequency  and  damping  ratio, 

as  given  by  (7.  5)  and  (7.  8),  are  plotted  as  fxmctions  of  the  ratio  of  depth 

to  radius  in  Figure  2  and  are  seen  to  approximate  unity  rather  closely 

when  the  depth  exceeds  the  radius  and  to  differ  therefrom  by  a  negligible 

amount  when  the  depth  exceeds  the  diameter. 

We  consider  as  a  niimerical  example  a  tank  10  feet  in  diameter 

and  of  equal  or  greater  depth.  The  period  given  by  (7.  5b)  is  1. 81  seconds, 

while  the  damping  ratio  given  by  (7.  7)  for  a  liquid  (e.g. ,  water  at  20°C) 

-2  -4 

having  a  kinematic  viscosity  of  10  stokes  is  4.  6  x  10  ;  the  latter 

figure  is  equivalent  to  a  time  of  7  minutes  to  damp  to  one-half  amplitude. 


•y  =  1.  12  1  +  (1  -— )  csch  (3.  68h/a) 
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The  displacement  thickness  given  by  (4*  5)  is  0.054  cm,  and  the 

corresponding  minimxnn  amplitude  necessary  to  render  the  laminar 

flow  tanstable  is,  according  to  (4.  7),  27  cm  (10.  6  inches).  We  note 

that  this  must  be  interpreted  as  a  (spacewise)  mean  amplitude  and  that 

it  is  based  on  the  minimum  transition  Reynolds  number. 

The  perturbation  pressure  on  a  vertical  elemen  f  the  wall  of 

a  circular  tank  may  be  calculated  from  {6.4b).  The  corresponding 

transverse  force  on  the  tank  is  given  by 

n 

,  Zir 

F  =  pgk  tanh  (kh)  C  cos  0  •  adO  (7.9a) 

J  o 

=  pgira^  I  ^  [  (ka)  ^  tanh  (kh)  cos  (cot)  {7.9b) 

where]  4  [denotes  the  maximum  amplitude  at  the  wall,  and  the  phase 
angle  has  been  arbitrarily  prescribed  as  zero.  Assuming  the  dominant 
mode  of  (7.4),  (7.9b)  yields  a  maximum  force  of  2700  lb  (oscillatory) 
for  a  sloshing  of  one-foot  amplitude  in  a  10-foot  tank  of  water.  This 
force,  according  to  (6.  6),  would  act  approximately  2.  7  feet  below  the 
mean  surface. 
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